Finite size effect on Bose-Einstein condensation 
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We show various aspects of finite size effects on Bose-Einstein condensation(BEC). In the first 
section we introduce very briefly the EEC of harmonicaUy trapped ideal Bose gas. In the second 
section we theoretically argued that Bose-Einstein(B-E) statistics needs a correction for finite system 
at ultralow temperatures. As a corrected statistics we introduced a Tsallis type of generalized B-E 
statistics. The condensate fraction calculated with this generalized B-E statistics, is satisfied well 
with the experimental result. In the third section we show how to apply the scaling theory in 
an inhomogeneous system like harmonically trapped Bose condensate at finite temperatures. We 
, calculate the temperature dependence of the critical number of particles by a scaling theory within 

O ' the Hartree-Fock approximation and find that there is a dramatic increase in the critical number 

of particles as the condensation point is approached. Our results support the experimental result 
which was obtained well below the condensation temperature. In the fourth section we concentrate 
on the thermodynamic Casimir force on the Bose-Einstein condensate. We explored the temperature 
"j^ ' dependence of the Casimir force. 
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O ' Bose-Einstein condensation (BEC) is a topic of high experimental [l|, 0, H, 0, d, 0| and theoretical [1, [l^, [ll|, 

interest. Within the last ten years several thousand works were done on this topic. Below a finite temperature a 
macroscopic number of Bose-particles come down to the single particle ground state. This phenomenon is called 
^ , Bose-Einstein condensation. Generally inter particle interaction is responsible for a phase transition. But the BEC 
' type of phase transition occurs entirely due to Bose-Einstein (BE) statistics. It is a consequence quantum statistical 
C^l [ effects of Bose particles. The theory of BEC was given in 1924 by Satyendra Nath Bose and Albert Einstein . But 
it took 70 years for the experimentalists to verify the theory P, [1, It is a confirmatory test of many body physics. 
Since 1995 a lot of BEC related experiments are being performed. The experiments opened up a lot of aspects of 
ultralow temperature physics and threw a lot of challenges to the theoreticians. 

Since the BEC type of phase transition does not occur due to the inter particle interaction rather it is a consequence 
"j'-y , of quantum statistical effect, the average separation(r) of the particles is required to be much larger than the s-wave 
■ scattering length(as) to probe this quantum statistical effect. Another necessary condition to probe this quantum 
^ \ statistical effect is that I is comparable to the thermal de Broglie wave length(AT). The P — T diagram of a Bose con- 
I ■ densate fall within the solid phase Q. So, the equilibrium temperatures and pressures in which the Bose condensation 
' O [ is achieved, correspond to the solid phase. Hence the Bose condensate is a different phase other than solid, liquid and 
^ ■ gas. It is a metastable state and in equilibrium it becomes a solid. Hence, to achieve Bose condensate one should 
take a very dilute Bose gas, otherwise the atoms will collide to form a solid. In the very low enough temperature the 
collision effects are less dominant. So one should go to the very low temperature to achieve the BEC. However, even 
at the low enough temperature, there is a possibility of two body collision. Due to the two body collision no energy is 
released out side the two body. Hence there is no possibility to form a molecule due to this type of collisions. In the 
5J] ■ low enough temperature there is finite possibility of three body collision. Due to this type of collisions the two body 
can release their energy to a third atom and they can form a molecule and eventually the whole condensate becomes 
a solid. Even if there is a wall the two atom can release their energy to the wall and the whole system can become a 
solid. For this reason the atoms are being trapped magnetically . For the trapping, the alkali atoms are good choice for 
their non vanishing magnetic dipole moment. The experimental studies of BEC started from 1970. The first studies of 
BEC were focused on the hydrogen atom which was considered because of its light mass, the most natural candidate 
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for realizing BEC. During the experiments hydrogen atoms were first cooled in a dilution refrigerator, then trapped by 
a magnetic field and farther cooling by evaporation took it close to the BEC. In the 1980s laser-based technique, such 
as laser cooling and magneto optical trapping were developed to cool and trap the neutral atoms. Alkali atoms are 
well studied to laser-based methods for their optical transition can be excited by available lasers and for they have a 
favorable internal energy-level structure to cool to very low temperatures. Once they are trapped their temperatures 
can be lowered further by evaporative cooling. Following the different cooling techniques the experimental teams of 
Cornell and Wiemam at Jila and of Ketterle at MIT succeeded in 1995 in reaching the very low temperatures and 
the densities to observe Bose-Einstein condensation in ^^Rb and ^■^Na respectively [i|,J2l- In the same year the BEC 
of ^Li was also achieved by Hulet group (3. Successively the BEC of ^He Q, ''^K [3 and ^"^Ci was 

achieved. The BEC of ^^Rb atoms which interacts attractively was also obtained by Jila group in 2001 lg|. In the 
last few years the experimental and theoretical studies of BEC-BCS crossover draw a lot of attention [171. Recently 
the experimental observation of strong quantum depletion in BEC has been achieved by MIT group [l8|. Since 1995 
the experimental study of BEC is unraveling a lot of new physics in the ultralow temperatures and is throwing a lot 
of challenges to the theoreticians. 

The theory of statistical mechanics as well as the theory of Bose-Einstein statistics and hence the theory of Bose- 
Einstein condensation is standing on the criteria of thermodynamic limit of the system. In the thermodynamic limit 
the volume(T^) of the system goes to infinity and the number(Af) of particles of the system goes to infinity such that 
^ is finite. However, according to the experimental setup of the BEC the system size is ^ mm^ and the number of 
particles is ^ 10^ P, 0, lii, 4, ^, Ji, lH ■ So, the criteria of thermodynamic limit is not satisfied. Hence we are interested 
in exploring the finite size corrections to the Bose-Einstein condensation. 



A. Theory of Bose-Einstein condensation for trapped ideal Bose gas 

There are two types of particles in the nature: fermions and bosons. Fermions have half integral spins and the 
bosons have integral spins. The Bose-Einstein condensation(BEC) has been achieved for ^^Rb, ^'^Na, ^Li, ''He, ^H, 
^^Cr etc. These are not the fundamental bosons like photons, gluons, Z-bosons etc. These alkali atoms are made up 
of even number of fermions (electrons, protons and neutrons). The spins of these fermions add up to make the atoms 
composite bosons. BEC is a quantum statistical effect of the bosons. It should be mentioned that the BEC of the 
photons is not possible because the number of photons is not a conserved quantity. Due to quantum statistics two or 
more identical bosons prefer to occupy a single quantum state. Below a certain temperature called the condensation 
temperature(To) a macroscopic number of bosons prefer to occupy the single particle ground state. A condensate is 
formed with the particles in the single ground state. For this reason the condensate gets a long range order and the 
condensate behaves like a macroscopic object. 

BEC has been experimentally achieved for harmonically trapped Bose gas. For the trapped system the condensation 
is achieved in the position space as well as in the momentum space. Let us consider a 3 dimensional system of 
harmonically trapped Bose gas, where all the particles are simple harmonic oscillators. Here all the particle are 
oscillating along a fixed point. Generally the oscillation frequencies along the three perpendicular directions are 
different. But, for the simplicity of the calculations we consider an isotropic case, where all the frequencies are the 
same. For the ideal case there are no interparticle interactions between the particles and the statistical mechanical 
behavior of the system can be well understood by a single particle hamiltonian. 

The energy eigenvalues of the single particle hamiltonian is 

3 

en^{n+-)hw (1) 

where n = 0, 1,2, 3.... and lo is the angular frequency of the oscillators. 

For the three dimensional harmonic oscillator the degeneracy of the state which has the energy {n + ^)fiLo is 
(n^ + 3n + 2)/2. The first term of the degeneracy factor is the bulk term and second term is the surface term. In the 
thermodynamic limit only the bulk term dominates. In the following section we shall show that the surface term will 
give the finite size correction. 
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B. Calculation of the condensation temperature(To) and condensation fraction 

We consider that the Bose gas is in equilibrium with its surroundings at temperature T. At and below Tq the 
chemical potential goes to {3/2)tiu). The total number of particles below Tg is obtained as 

= r fill . dn = [|^]\(3) (2) 
At To there is no particles in the ground state. So, we must have 

N = [^]^C(3) (3) 
Prom the above equation we get the condensation temperature as 

'^" = tW 

Let us now define the thermodynamic limit of the harmonically trapped Bose gas. In the thermodynamic limit oj ^ 
and N ^ 00 such that Nu^ is a finite constant. From the eqn.(2) and (3) we get the expression of condensed fraction 
as 

f = l-f = 1-[T/TJ3 (5) 

From the above equation wc see a macroscopic occupation of the particles to the lowest single particle energy state 
is possible below a finite temperature T^. This is a consequence of Bose-Einstein statistics and this phenomenon is 
called Bose-Einstein condensation. 



C. Specific heat 



The total number of particles above Tg is obtained as 

,2 , 



Prom this relation we get 



The total energy of the gas is obtained as 

,2 , 



/■OO „2 /O h.rp 

E= / nhw , ' — -dn = 2,kT[—Ygi{z) 



for T > To- From the above eqn.(8)and (7) we get the specific heat as 



(6) 



(7) 



(8) 



kT kT 

a = m—rg4{z) - 9k[—]'gUz)/92{z) (9) 

At and below To we have z — 1. So the specific heat for < T < To is 

kT 

a = 12k[^fg4{l) (10) 
From the eqn.(9) and (10) we see that there is a discontinuity in C„ at T = Tq. The amount of discontinuity is 

AC, = -Q^^Nk (11) 

This discontinuity (or the discontinuity in the derivative) in specific heat signalize that BEC is a phenomenon of phase 
transition. 
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D. Shifts of To 

In the previous section we stated that the degeneracy factor of a state of the 3-d isotropic harmonic osciUator with 
energy (n + 3/2)huj is (n^ + 3n + 2). In the thermodynamic limit only the first term of the degeneracy contributes. In 
the following calculations we will see that the second term on the degeneracy factor will give a finite size correction. 
The corrections due to the third term on the degeneracy factor and that due to summation integral conversion are 
negligible. Now including the second term of the degeneracy factor we get total number of the excited particles below 
the condensation temperature as 

/•°° 1 f°°3n 1 

Jo 2 e"'*'^/'^'^ -I Jo 2 e"'^'^/'^^ - 1 

Due to the inclusion of the second term of the degeneracy factor we get the new condensation temperature Tc such 
that 

Comparing the eqn.(3) with the eqn.(13) we get a shift of the condensation temperature Tc — To, such that 

at; Tc-Tq C(2) 1/3 ^ Q 73Ar-i/3 ^4) 

To n 2[C(3)]2/3'' 

From the eqn.(14) we see that the finite size induce a lowering of the condensation temperature. 



II. MORE ACCURATE THEORY FOR BOSE-EINSTEIN CONDENSATION FRACTION 

More accurate theory for Bosc-Einstcin condensation fraction was discussed in [20j . 
In the thermodynamic limit we have the B-E statistics as 

The mass of each particle be m. The length and the volume of the system are L and V{— L^) respectively. Thermal 
de-Broglie wave length of a single particle is \t ~ \J ^7^, where k is the Boltzmann constant. Average separation of 
particles is r= (^)^/"^. In the classical limit we must have 



I 



. »1. (16) 
At 

From eqn.(16) we can easily write 

kT 

(17) 

At sufficiently low temperatures when ^ 1, the gas becomes degenerate and quantum correction is necessary. BEC 
occur at the onset of this degeneracy. So, for T ^ To we have ~ N'^/'^. However the condition of thermodynamic 

limit(i/AT 00) is also valid at this temperature. For finite system the thermodynamic limit is realized as L/Xt 3> 1. 
However, at sufficiently ultralow temperatures L/\t can be comparable to 1. In this situation, the B-E statistics 
needs a correction. 

The usual condition of thermodynamic limit is not properly satisfied in the case of the experimental setup of BEC 
of 3-d harmonically trapped Bose gas[ll, H, [EJ3, S S [3|- the semiclassical approximation the number density at 
the distance(r) from the center of the trap isQ 
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where At — y ^7^; and 53(2;) — x + jl?!"^ + + ... is Bose-Einstein function of a real variable x. n{r) 



,2^2 



e 2fcT ^ where -^muj r is the trap potential and uj is the angular trap frequency. The length scale of this 3-d 



trapped Bose gas is Lt ~ \J Putting this expression of length in eqn.(17), we get the conditions of classical 

limit as ^ N^^^, and the thermodynamic limit(LT/AT 00) as ^ 00. Eventually we get the condensation 
temperature for this system as To ~ ^N^/^. More precisely the thermodynamic limit for this system is realized as 
to — > 0,7V 00 and Noj'^ = constant. From the above semiclassical expression of number density we have @ 



With the consideration of thermodynamic limit, total number of excited particles for T < To are Nt = [■^]'^C(3) 
and the number of condensed particles are No = N[l — (^)"^]. In the experimental setup the typical value of[3| TV, 
u! and m are of the order of 50000, 2645s~^ and 52 amu. With these experimental parameters, at To the length 
scale as expressed in eqn.(19) is Lt ~ 56 x 10~^mm. and Xt„ = 2.88 x 10~'^mm. and their ratio is Lt/Xt ~ 32.2. 
However, for T ~ a large fraction of particle come down to the ground state and the length scale of the system 

becomes ~ \J ~ 6.76 x lO^Snm. At these ultralow low temperatures, the thermal de Broglie wavelength becomes 

comparable to the system size. At these ultralow temperatures the usual theory of statistical mechanics of finite 
system is not properly valid. So, for ^ ^ T < To, we seek an ultralow temperature as well as finite size correction 
to the B-E statistics. 



A. Correction of B-E statistics from qualitative point of view 

To quantify the correction arising from ultralow temperatures let us start from Tsallis statistics f2l'|. The relative 
probability that E be the total energy of a system is given by Boltzmann factor e~^/^^ . In the Tsallis statistics 
Boltzmann factor is replaced by (2+(g-i)£;/A;T)i/(g-i) ' where g is a hidden variable. It is easy to check that as 1, 
we get back Boltzmann factor. However, for the Bose particles we should start from the Tsallis type of generalized 
Bose-Einstein statistics which is expressed as|22i] 

(20) 



Once again it is easy to check that, in the generalized B-E statistics, as q 1, we get back B-E statistics. 
However, Tsallis statistics^] is applied to equilibrium 2^^ [2^ . [2^ as well as to nonea uilibrium [26l. [27l |28| systems. 
Tsallis statistics is re-derived as one of the superstatistics[23| of a nonequilibrium system. In this theory the hidden 
variable(g) of Tsallis statistics is equated with system parameter and q is no longer a variable. In the theory of 
superstatisticsjl^ 2/(q — 1) is redefined as effective no. of degrees of freedom. In the theory of dynamical foundation 
of nonextensive statistical mechanics, this effective no. of degrees of freedom is equated as[28| {2i — q)/{q—l). However, 
for finite equilibrium system we equate q with a system parameter. We will equate g — 1 with Xt/ Lt so that in the 
thermodynamic limit(LT/AT ^ 1) of finite system we can go to the usual Bose-Einstein statistics. FIG.l shows 
experimental and theoretical plot of condensation fraction with temperature. In this figure we see that a negative 
shift of condensation fraction is necessary to satisfy the theoretical plot with the experimental data. In the theoretical 
plot there is finite size correction and the correction due to two body interactions. We shall see that the generalized 
B-E statistics with (g — 1) cx |y will give rise to a significant negative shift of condensation fraction with temperature. 
This significant shift along with finite size correction and the correction due to the two body interaction might satisfy 
the experimental results. For the trapped Bose gas we shall start from eqn. (20) and put g - 1 = to get 

^ (21) 



il + ^-fcf^J 



where a is an arbitrary constant. This a is to be determined from the experimental result. 
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FIG. 1: Condensation fraction (■^) versus temperature(-p) plot. The thick hne follows from equation(26). The dashed line 
corresponds to the thermodynamic limit and excludes all the correction terms in equation(26). The dotted line corresponds 
to the finite size correction and the correction due to interaction excluding the ultralow temperature correction term in equa- 
tion (26). All the theoretical curves are drawn according to the following experimental parameters. The dotted points are 
experimental points of Bose-Einstein condensation of ^^Cr whereQl N = 50000,LJ = 2645s"\ m = 52a .m.u.,a = 105a_g and 
To ~ imnK 



B. Corrected and generalized B-E Statistics 



We are considering a finite system. For tiiis finite system the condensation temperature would not be To- For this 
system the condensation temperature is T^ which must be close to Tq. For a system of Bose gas in isotropic harmonic 
trap, the single particle energy levels are Cj = (| + j)^, {j = 0, 1, 2, 3, ...). For T < Tc, we have /z = f^"^- So for 
T <Tc, from eqn.(21) we get the average number of particles in jth state as 

^ (22) 



3 lat _ 1 



where t = j^. From the above eqn.(22) we get the expression of the total number of excited particles below the 
condensation point of a 3-D Bose gas isotropic harmonic oscillators as |20| 



3f2 f,f2 

Nt = t^C(3) + — C(2) + — C(4) (23) 
2 a 

The gas being very dilute there should be a correction term in the expression of Nt due to two body scattering. 
This correction within Hartree-Fock(H-F) approximation has been discussed in [s^. According to H-F approximation, 

the correction term to the above Nt is 3 x 1.326 , ° nI/^ — 4.932 . ° t^/^, where a is the s-wave scattering 

wh/rriLj wh/rnuj 

length. Introducing this term in eqn.(23) we get the more corrected expression of number of excited particle as 

Nt = t^Ci^) + ^C(2) + 4.932-^i^/2 + —((4). (24) 
2 ha 



At T — Tc, all the particles will be in the excited statesjl, Q. So at T = Tc or at tc — the number of excited 
particles is equal to the total number of particles. So, we must have 

Q+2 af2 

N = tlCiS) + ^C(2) + 4.932^=i^/2 ^ ^^(4) (25) 
2 h a 
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FIG. 2: Plot of r'{t) with Temperature(t) in units of This plot is followed from eqn.(29). The doted line shows the 

theoretical value of the ratio r'{t) in thermodynamic limit when uj —> 0. 



Corrected and more accurate condensation fraction 



In the thermodynamic limit of the ideal trapped Bose gas[l, Nt = ^'^C(3) and its condensation temperature To 

is such that[l, [§] to = , where t^, = Comparing this expression of to and tc of eqn.(25) we see that 

tc < to and there is a shift St^ — — to of condensation temperature due to the inclusion of finite size correction 
term, correction term due to two body scattering and due to the ultralow temperature correction of B-E statistics. 
For t < tc, from eqn.(24) we get the fraction of number of particles in the ground state as 



No 
N 



T 



N 



C(2) +4.932 



f/2 



6t' 



(26) 



D. More accurate To shift 



However, in the thermodynamic limit with no inter particle interaction, the expression of this fraction would be 

Vol _ 1 

w\t-l - i - 
temperature as 



^]t-l = 1 ~ ['t]'^- N^ow, from eqn.(26), due to the correction terms we get the fractional change in condensation 



^ ^ _ 1 ^9P, " /yi/e ?C(4)_ 1/3 ^27) 

To 2[C(3)]V3^^ '-'''yftA^'' a[C(3)]^/3 ^ 



From the first term of the eqn.(27), we get the Tc shift due to the finite size correction as [l£ 



stl 



to 



~ 2[a3)]V3 ^~^^^ = -.7287V"i/3 and for 50000 particles we get = -1.97%. From the second term of the eqn.(27), 

we get the Tf, shift due to the correction of two body interaction as [1, [^, [131 % — — —1.326^=2 — iV^/^ = —6.61% 
for uj = 2645s~\ a = 105as and for N=50000 Q. For this setup to = [^1^^^ = 34.65. From the third term of the 
eqn.(27) we get the Tc shift due to ultralow temperature correction of B-E statistics as = — ^^^^'3^^2/3 N^^^"^ = 
— ^5.19% for 50000 particles. According to the experiment should be 10%. This 10% tc shift is achieved if we put 
a = 1.48 in eqn.(27). 
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E. Comparison between the specific heats calculated by the generalized B-E statistics and B-E statistics for 

the Trapped Ideal Bose Gas 

Well below Tc, as t approaches to 1, the B-E statistics needs more corrections. Below Tc, as a comparative study 
of the two statistics, let us calculate the ratio of specific heat from the two statistics. As a comparative study we can 
take ideal gas. The total energy of the system would be 

E(t) - (huj) / 4r ■ dj 

Jo 2 [1 + 



As i > 1, E{t) -> 3SwC(4)i'* as expected from the B-E statistics d Q. 

Let us denote the ratio of specific heat(C„(t) = ■^^E{t)) calculated from the corrected B-E statistics and from 
the B-E statistics as 

Obviously as t ^ oo, r' 1. With a = 1.48 numerical plot of r'{t) at very low temperatures (1 ^ i < tc) is shown in 
the FIG. 2. 

In FIG. 2 we see that at high temperature t ^ 1, the specific heat behaves well according to our familiar law. In 
this figure we also see that specific heat for this finite system becomes negative at i < 5. In this range of temperatures 
Lt ^ At and the theory of statistical mechanics is not valid. In spite of that, due to nonextensivity the appearance 
of negative specific at these temperatures is not surprising theoreticallv[3l| and experimentally [s^ . 



F. Conclusion 



That B-E statistics needs a correction for finite system at ultralow temperatures is justified theoretically. But we did 
not theoretically justify why Tsallis type of generalized B-E is necessary for finite system at ultralow temperatures. 
However, Tsallis type of generalized B-E statistics with our redefined parameter satisfy experimental result. We 
argued that the correction over the B-E statistics is 0{\t / Lt)- Smaller the ratio smaller is the correction. However, 
this ultralow temperature correction is valid only for a finite system. How a is to be determined theoretically remains 
an open question. 



III. FINITE TEMPERATURE SCALING THEORY FOR BOSE-EINSTEIN CONDENSATE 

Finite temperature scaling theory for Bose-Einstein condensate was well discussed in [33] . For the ultracold gas 
the interaction is characterized by s-wave scattering length a^. Atomic interaction as well as the value of scattering 
length(as) can be controlled in Feshbach resonance [s^l- Stability and collapse of Bose gas with negative scattering 



length has been observed in the clouds of ultracold ^Li[35[ and *^Rb[36|, l37[. If the interaction is attractive (a^ < 0), 



the gas tends to increase the density of the central region of the trap. This tendency is opposed by the zero-point 
energy and thermal energy of the atoms. If the number of atoms is greater than a critical number(A^c)i the central 
density increases strongly and the zero-point and thermal energy are no longer able to avoid the collapse of the gas. 
As the temperature (T) goes to zero, all the particles come down to the ground state and the system is well described 



by the ground state wave function 5*0 (r) = y j3~3j2^ ^ iii the position(r) space, where I = ^hjrnu} is the length 
scale of the oscillators. At T = 0, the density of the condensed particles is described as no(r) =| ^'o(r) P= ;3^/2 e~"^ . 

^^2^2 y 

In absence of collision the number density of the excited particles isiS] UTij) = -^g^{e~ '^^bt )^ where \t = \j ^ , 

and (73 [x) — X + x"^ 12^/"^ + x^ /i^/"^ -I- ... is the Bose-Einstein function of a real variable x. 

Now, if we allow two particle attractive interaction, the particles come closer and the length scale of the system 
would reduce by a factor v (< 1) such that the ground state wave function would be 'I'o('') — ^ ^2;3°3/ 7e~'' ^^"^ ' and 

the thermal density of particles would be ^^(r) — ^s^i gi (e ^"'^''bt ). 
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Let us consider the interaction potential as ymt(r) = g5^{r), where g = — ^'^^ ° is the coupling constant and 
a = —Qs is the absolute value of the s-wave scattering length. For a dilute gas we must have j ^ 1. The typical 
two body interaction energy for N number of particles is ~ N'^g/21^. For this interaction, the gas tends to increase 
the density of the central region of the trap. Well below the BEC temperature (Tc), i-e. for T 0, this tendency is 
resisted by the zero-point energy (~ Nhw) of the atoms. In the critical situation the typical oscillator energy must 
be comparable to the typical interaction energy. So, at T = 0, we must have Nchuj ^ N^g/2l^. From this relation we 
can write ^ 1. 

For < T < Tc, the typical total energy of the system is N^/^fno. So at these temperatures we must have 
Nc'^hw ~ N^g/2l^. From this relation we have ^ - [^]^^^ > 1. 

However, near Tc the length scale of the system is@| Lt^ ^ ^\/^f^ ^ IN^/^ . So, near the condensation temperature 

we must have Nc^^huj ^ Ncg/2L^^. From this relation we have ~ [^]^ » 1. 

In the following subsection we shall explicitly calculate by a scaling theory within the Hartree-Fock(H-F) 
approximation. We shall explicitly show the temperature dependence of the critical number(7Vc). 



A. Energy expression of the Bose particles within Hartree-Fock approximation 

Although the problem of an attracting Bose gas was discussed by many authors [olIToj. yet the temperature depen- 
dence of the critical number of particles has not been explicitly explored. H-F approximation for Bose gas has been 
discussed in Q. For this interacting Bose gas, the average occupation number fii as well as the single particle wave 
functions 0i are determined^ by minimizing the grand potential. Since the entropy and the total number of particles 
depend on the occupation number n^, the single particle wave functions {(f>i} are simply obtained by minimizing the 
energy functional with proper normalization constraint j <Py \ (jji \'^— 1 for each i,8]. Within the H-F approximation 
we have the expression of energy functional asfs'l 




-KF(r)no(r) + V{Y)nT{v) + ^-nl{v) 

+2gno{T)nT{Y) + gnl{T)] (30) 



B. scaling argument 

In our scaling theory we introduce a variational parameter v which would fix the width of the system as well as would 
take a proper choice of no(r) and ^^(r). So, in this method, the choice of is denoted by the variational parameter 
u. The minimum of the energy functional for a certain choice of would corresponds to the equilibrium of the 
system. So, for equilibrium condition of the system, the energy functional would be minimized with respect to the 
variational parameter (z^). In FIG. 3 we shall see that above a critical number of particles(iVc), the energy functional 
will have no stable minimum. If the energy functional as well as the grand potential has no stable minimum, the 
grand potential as well as the energy functional would arbitrarily decrease until the width(i^) of the system becomes 
zero. Under this condition the system is said to be collapsed. 

The scaled form of ^^(r) is obtained from the variation of B-E statistics such that n(p, r) — — 7-2 , 

where p is the momentum of a single particle. With this variational form of statistics we have the total number of 
excited particles as iV = ( '^n^}" )^C(3)- For T < the total number of particles in the ground state is uq = N — Nt = 
N{1 — ^)^. From the variation in statistics we evaluate the energy functional of eqn.(30) in units of Nhuj as[3^ 

X,M=ci(i)(4 + ^^)-^ (31) 

where t = T/T^. c,{t) = |(1 - t^) + I^Q}^ and c,{t) = ^^^[1 - t^]^ + ^^^i£^t^n[l t^] + 
such that S' = Erj=i (.j)3/2(',+,-)./2 ~ 0.6534. Putting f = 0.0066 and i = 0.8, in eqn.(31) 
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FIG. 3: Plot of energy per particle (Xo.8(i^))in units hu) with the variational parameter u. For the thick line, — 27.78. For 



the dashed line, 



20. For the dotted Une, 



12. All the lines follow from eqn.(31). 



we get the FIG. 3 for various number of particles. In FIG. 3 we see that Xt(v) has a stable minimum and an unsta- 
ble maximum below a critical number of particles. At the critical number of particles the minimum and maximum 
coincides. Above the critical number of particles there is no stable minimum and its energy arbitrarily decreases to 
— oo and its width becomes zero which signalize the collapse of the gas. 



C. Evaluation of the critical number 

For the critical number of particles we must have — ^tt^ 0. From these critical conditions we get the 



expression of for i = as 



_2 fn-r 7^ — n .jc 



-^ = ^[^,-^l]=Omi (32) 

The same result was also obtained in 'so']. However, the experimental value is 0.459 [l^. From the same critical 
conditions we get the expression of for < T < as 



x(l-t3) + 1.636t^/^][^]V4_!__ (33) 



Eqn.(33) as well as FIG. 4 represent the relation between the critical number of particles and the temperature. 



However, from this figure and according to eqn.(33) we get — 27.89 for t — 0.8. From the FIG. 3 we see that 



N, 



; - 27.78 for t = 0.8. 
The same critical conditions give the expression of for T = Tc as 



(34) 
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FIG. 4: Plot of the critical number Nc in units of | with temperature(r) in units of (Tc). Here ^ = 0.0066 [3. This plot 
follows form eqn.(33). 

Comparing the eqn.(33) and eqn.(34), we see that there is a dramatic increase of as the condensation point is 
approached. It would be interesting for the experimentalists to do the experiment near Tc and to verify this dramatic 
increase of the exponents of a/ 1 in the expression of Nc- 



In the FIG. 4, we see that the temperature dependence of Nc is not significant for < T < O.STc. If the experiments 
of the collapse is performed within this range of temperature, the zero temperature theory should well satisfy the 
experimental results. Here we have assumed that the radius of the thermal cloud changes by the same factor as that 
of the condensate cloud due to the fact that Lt is proportional to I. Initially we qualitatively explained the physics 
of the collapse of the attractive atomic Bose gas. Then we semi qualitatively estimated Nca/l for the various range 
of temperatures (0 < T < Tc). Finally we calculated Nca/l by a scaling theory (variational method) within the H-F 
approximation. Our calculation supports the semi qualitative estimation of Nca/l. It proves that exponent of a/ 1 in 
the expressions of Nc for the different range of temperature does not depend on the scaling assumption. The same 
temperature dependence of Nc as well as the same type of exponents can be obtained by a quantum kinetic theory. 

IV. CASIMIR EFFECT ON BOSE-EINSTEIN CONDENSATE CONFINED BETWEEN TWO SLABS 

Casimir effects on Bose-Einstein condensate confined between two slabs were well discussed in ji^]. Vacuum 
fluctuation of electromagnetic field would cause an attractive force between two closely spaced parallel conducting 
plates. This phenomenon is called Casimir effect and this force is called Casimir force [4l|, |42, . Before going into 
details we first estimate the Casimir force intuitively. Since this force causes due the vacuum(quantum) fluctuations 
of the electromagnetic field there will be a term h in the expression of vacuum energy. As we are considering the 
vacuum of electromagnetic field, there must be a term c(velocity of light) in the expression of vacuum energy. Since 
the quantized energy of the electromagnetic field is h times frequency we can conveniently write the expression of 
vacuum energy as he/ L, where L is the plate separation. Since the fluctuation are confined between the plate, the 
Casimir energy is expected to be proportional to the area(A) of the plate. For the proper dimension, the Casimir 
energy would be ^ Ahc/L^. So the Casimir energy would be ~ Ahc/L^. This is an estimation of the Casimir force. 
In the following we show an exact calculation of the Casimir force. 



D. Conclusion 
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Let the plates be kept parallel to the x-y plane. They are separated L along the z- axis. The vacuum energy of an 
electromagnetic field \s E = Sk^^iWk, where WkS are the angular frequencies of the photons. The vacuum energy, in 
term of the wave number(A:2:, ky,kz) can be written as 



EiL) = 2Ax^^hcE^^^J^ ^kl + ^27:k^dk^ (35) 

where k\ — k^ + ky and kz — rm/ L. The factor 2 of the above equation comes due to the fact that photon has two 
polarizations. With the substitution k\ + ^^^jl- — and A as the cutoff of t we can recast the eqn.(35) as 

From the definition of Casimir force (Fc{L) — —-^[E{L) — E(oo)]) we get the expression of Casimir force as 

Ahcn^ 3_ Ahcn^ A A 

W)--^^S„^o" - -^I^^(-3) - - 2^ (37) 

where C(~3) — is obtained from the analytic continuation technique. This expression of the Casimir force is valid 
only at zero temperature. This force has been measured experimentally [i^. 14^. 

However, Casimir effect can be generalized [i^ for any range of temperature and foi^any dielectric substance 
between two dielectric plates. It has also been generalized for thermodynamical systems 43]. Casimir force for this 
kind of systems has recently been measured At finite temperature T, the definition of Casimir force is generalized 
as gil 

i^e(^,i) = -^[^^T(i)-^^T(oo)] (38) 

where VtT{L) is the grand potential of the system confined between two plates separated at a distance L. 

We consider the Casimir effect for a thermodynamical system of Bose gas between two infinite slabs. Geometry 
of the system on which some external boundary condition can be imposed is responsible for the Casimir effect. 
Thermalized photons (massless bosons) in between two conducting plates of area A at temperature T gives rise to 
the Casimir pressure [49| 

FAL) TT^hc 16(fcr)4L4 ^fj^ 

~ 240L^[^ + 4(4^] ^^'^fcTL^l 



where k is the Boltzmann constant, c is the velocity of light and L is the separation of the parallel plates. At T ^ 0, 
Casimir pressure 1 
high temperature 
independent of fi. 



Casimir pressure becomes — ^aql^^ ^"^^ only the vacuum fluctuation which contributes to the Casimir pressure. At 
high temperature i.e. for 0, the Casimir force for photon gas goes as and has a purely classical expression 



A. Calculation of the B-E condensation temperature 

Let us consider a Bose-gas confined between two infinitely large square shaped hard plates of area A. The plates are 
along x-y plane and they are separated along z- axis by a distance L. For the slab geometry, \J~A ^ L. We consider 
that our system is in thermodynamic equilibrium with its surroundings at temperature T . At this temperature the 

thermal de Broglie wavelength of a single particle of mass m is A = \^ -imkT • ^'^^ thermodynamic limit, ^ ^ 1. For 

this system the single particle energy is ^{px,Py,j) = ^ + |^ + '^2mL'^ ' ^here px and py are the momentum along 
X-axis and y-axix respectively and j — 1, 2, 3, However in the thermodynamic limit the single particle energy can 

2 2 2 

be written as e{px,Py,Pz) — §^ + + where Pz is the momentum along z-axis. 
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Considering the thermodynamic hmit the total number of thermally excited particles can be written as 



OO />00 /"OO 



e kT — 1 

where is the chemical potential and V is the volume of the system. For this system the BEC temperature can be 
calculated as 

B. Finite size correction to the grand potential 

Let us now introduce the finite size correction. The ground state energy of our system is [g — 2rnL-i \ ■ The average 
no. of particles with energy ep p j is given by — 5 5 ^ where /i' = (/^ — g) < for bosons. At 

and below the condensate temperature — > 0. For this bosonic system we have the grand potential 



a{A,L,T,^i') 



g kT g 2mkT g 2mkT g 



frp-»('^(T)'[/'+2/]) 



(42) 



where A = 1 / 2^kT ■ Integrating over and pj, we get 



L,T, = -^^^ E E — - 2/-^ + ~ 3^ + (43) 

i=lj'=0 

Since <C 1, higher order terms of the above series would not contribute significantly. From Euler-Maclaurin 
summation formula we convert the summation over j' to integration. The Euler-Maclaurin summation formula for a 
smooth function f{k) is denoted as 

- £ + - T2^'W + - 30^0^ ■■■■ 

Now from eqn.(43) and eqn.(44) we have [i^ 

17(Ai,T,Ai') 

i— 1 

f^^/^^ + o([^]^)] 

A{kTYm.L . . 1 , , ttA , 
= -^^[2A^t(") - 2^^(") + 2L^i(")] 



(45) 



where z = /'''-^ is the fugasity and gi{z) — z + ^ + ^ + .... is the Bose-Einstein function. From the above equation 
we get the total number of particles as 
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FIG. 5: Casimir force(Fc(r, i)) versus temperature(r) plot. Fc{T,L) is in units of and T is in units of To. Equation 

(47) corresponds to the range ^ < 1. Equation (48) corresponds to the range 1 < < 1.1. 



N 



AkTm L 1 ttA 
— 3 (z) Qt{z) H qi (z) 



(46) 



In the thermodynamic hmit of a system, as T < To, z —> 1. For a finite system this can not happen, otherwise 
the correction terms in the above expression would be infinite. Instead at T > To, z ~ 1. Taking only the first 
correction term in the eqn.(46) we have g\{z) = -~ln(\ — z) — \N' {T)g^{z) — A^C(3/2)] ^pfifrp^ — —InAz, where 

N'{T) = 4^[^C(3/2)] and Az = 1 - z is a small change in the fugasity at T > To. Now putting z = 1 in the 
expression of Az, we get Az = e"^^'^^-''/^)^/^^'^^)^), where AA^ = N'(T) - N. We see that in the thermodynamic 
limit(i oo) Az — and when L is finite such that L/ X ^ 1, we have z ~ 1 at T > Tq. 



Calculation of the Casimir force 



Let us now calculate the Casimir force. At T < To we put fj,' ^ or z ^ 1. Now, from the eqn.(45),(41) and (38) 
we get the expression of Casimir force for < T < Tc as [40| 



(47) 



Above the condensation temperature /i' < or z < 1. However, for T > To, z ~ 1. Putting T = To + AT in 
equation (45) and from the definition of Casimir force we have [40] 



F,{T,L)^-N 



(48) 



where < ^ 1. Hence, for T > To, the Casimir force weakly depends on temperature. 



D. Casimir Force for T > To 



Let us now calculate the Casimir force at T ^ To. At these temperatures z ^ 1. So we can approximately 



write (?i(z) !v z + From the first term of eqn.(46) we have (z) 



N2-iTh^2X 
AkTmL 



z + z^/(2-\/2). For this range of 
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temperatures we can write e^^/^'^ = ze '^^ k, z. For convenience, we replace /i' by /i and recast equation (43) as 

m 

o(^,L,r,rt.-j^i:j:— (^--+3^.^) (49) 

where we use the formula X]^^-oo ^ ~ X^riL-oo ^ ^'^^ T" ~^ expression of above Casiniir 

potential(the third term of the grand potential) we can put e'^/^'^ = z <^ 1 and can take i = 1 and j = 1 as the 
leading term to contribute in the Casimir potential. So, for T ^ To we have the Casimir force as [13] 

F,iT,L) = ^-^ e (50) 

Now we see that in the classical limit(r ^ To) the Casimir force vanishes as e"*^^. It is interesting to see that the 
Casimir force for T ^ To is long ranged(power law decay) and for T < To is short ranged(exponential law decay). 

The changes of Casimir force with temperature for the range < T < To and for the range T > To is shown in 
FIG. 5. 



E. Casimir force for trapped geometry 

Casimir effect on BEC for trapped geometry was discussed in (soj. We consider a system of ideal Bose gas which 
behaves as isotropic trapped harmonic oscillators along two perpendicular directions and along the other perpendicular 
direction it behaves like particles is 1-d box. At thermodynamic limit the system size must be greater greater than 

the thermal de Broglie wavelength ^ 2mkT particles. For a slab geometry of a 3-d system the length scale along 

one direction is much less than that along the other two perpendicular directions. Similar type of analysis for this 
system gives the expression of the Casimir force as [HO] 

Fo{T,L) = ~N[^]'/^^ forT<To (51) 
lo mL'^ 

Here the exponent 5/2 of T/Tq is a consequence of the trapped geometry. 



F. Conclusion 

That vacuum fluctuation causes Casimir force is well known [4l[. Critical fluctuations also causes Casimir force [43j. 
The Casimir force calculated here is neither due to vacuum fluctuations nor due to critical fluctuations. It is due 
to quantum fluctuations. These fluctuations are associated with the commutator algebra of position and momentum 
operator as well as with the commutator algebra of bosonic annihilation operator(di) and creation operator (dj) such 
that [cLi, a}-] — Sij, where i,j represent the single particle energy states. For T <To & macroscopic number of particles 
come down to the ground state. The quantum fluctuations die out due to the macroscopic occupation of particles in 
a single state. That is why the Casimir force dies out at T < To . For, T ^ To the Bose-Einstein statistics becomes 
classical Maxwell-Boltzmann statistics and thermal fluctuations dominates over the quantum fluctuations. For this 
reason the Casimir force dies out at T 3> To . For our system, the reduction of thermodynamic Casimir force with the 
power law will show the signature of the Bose-Einstein condensation. 



APPENDIX A: ZETA FUNCTIONS FOR NEGATIVE INTEGERS 

Zeta functions for the negative integers are obtained from the analytic continuation technique [sij . However, here 
we will present a simple technique. Our technique is as follows. 

Let us consider a function / of a variable t, such that /(t) = eTTT- We will expand this function in two different 
ways and will collect the coefficients of the powers of t. Equating the coefficients we will get the zeta functions for 
the negative integers. 
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fit) 



1. Expansion of f{t) 



t 



e* - 1 



= te' 






he-2* 


+ + 






i3 *4 

3! ^ 4! 

2^*3 




22*2 


2^*4 


~2r ' 

32*2 




~W 


2! 





= t [ 1-iH 
+ 1 -2t + 
+ 1 -3t + 

+ ] (Al) 

We see that the coefBcient of t'' of the above expansioii(Al) goes as (— 1)^+^ '^[k-iy} ^^^^ 

/(fc) = Sr., (-^)y5-^^ + ^^ *^ (A2) 
2. Another expansion of f{t) 



(- -|- 2! 3! ^ 4! ^ ••• ^ 2! ^ 3! ^ 4! ^ •" 

t t , 1 

+ 2! + 3!+4!+-) 
4! 3! 4! ^ 2! 4! 3! 4! 3! 4! 3! 4! 

1 ,1 1, , , 1 1 1 , 1 , , 1 1 1 11 1 , . , 

= ^ - 2^ + (4 - 6^^ + (-4! + ^2! 3! " ^2! ^ + 5! + W ^ ~ 'W^- ^ W^' ^ - ^^'^ 
The coefficients of ** the above expansion(A3) can be written in terms of BernouUi's number(i?;;) such that 

/(fc)=S-o|f<'- (A4) 
From the above expansion(A3) we see that Bq = 1, Bi = —1/2, B2 = 1/6, B3 = 0, i?4 = —1/30, B^ = etc. 

3. Relation between the Bernoulli's numbers and zeta functions 

Let us now equate the coefficients of t'^ of the two different expansions of f{t). Equating the coefficients we get 

Ci-k + 1) B, 

^ ^' (fc-1)! - fc! ^^^> 

Prom the above eqn.(A5) we get 



a-k) = (-1)'^^ (A6) 



Now wc sec that C(0) = 1 + 1 + 1 + 1 + .. = Bi/1 = -1/2, C(-l) = 1 + 2 + 3 + 4 + .. = -B2/2 = -1/12, C(-2) = 
12+22+32+42 + .. = 53/3 = 0, C(-3) = 1^+23+33+43 + .. = -B4/4 = 1/120, C(-4) = l''+24+3''+44+.. = B5/5 = 



17 



etc. Although our first expansion of f{t) is not mathematically sound yet the results obtained for the Ci^k) are the 
same as obtained from the analytic continuation technique fsij. 
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